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ABSTRACT 

Simultaneous measurements of waves, and vertical and 
horizontal water particle velocities were made at the 
breaker-line within the surf zone using a capacitance type 
penetrating wave staff, a pressure wave gauge, and an 
electromagnetic current meter. Wave' measurements were also 
made at seaward and shoreward locations. The wave energy- 
density spectral components were converted to velocity 
spectral components using linear wave theory. These 
computed values compared well qualitatively with the meas- 
ured velocity spectra. Quantitatively, the results showed 
that linear theory underpredicted wave-induced horizontal 
velocity spectral components by about 50 percent at the 
frequency of peak energy. The coherence values between 
waves and horizontal velocity were high, ranging above 0.75* 
The phase angle computation showed the calculated velocity 
components leading the measured velocity components by an 
average of 20 degrees, indicating an unstable wave crest 
leading the particle motion in the body of the wave. 
Probability density functions were computed and compared to 
Gaussian and Gram-Charlier distributions using the chi-square 
goodness-of-f it test. The Gram-Charlier distribution 
qualitatively gave the better fit to the data. 
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I. INTRODUCTION 



A. NATURE OP THE PROBLEM 

The problem is to describe the kinematics of breaking 
waves within the surf zone in terms of the water surface 
elevation. The theories which have been developed for deep 
water waves can be applied with some degree of certainty, 
and can be tested in a laboratory situation. These theories 
generally hold up to the time that the wave begins to shoal, 
and with somewhat lesser accuracy, throughout the shoaling 
process up to the point of breaking. At the breaker point, 
however, the prediction of wave kinematics becomes more 
difficult. The type of breaking w ave, plunging, surging, or 
spilling, depends upon deep water characteristics, bottom 
type, and bottom slope, hence requiring a different descrip- 
tion of the kinematics within the wave for each breaker 
type. Physical limitations preclude the construction of a 
realistic laboratory model of the surf zone. And finally, 
the hostile character of the surf zone complicates the 
taking of ill situ measurements to the degree that relatively 
little data exist from within the surf zone. 

The most forthright approach to the problem of describing 
the kinematics of breaking waves in the surf zone seems to 
be through direct measurements. Advances in instrument de- 
sign have led to simple, sturdy, measuring devices with 
rapid response time which can sense the small scale as well 
as the large scale features. By making the assumption that 
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the processes are stationary, this assumption being valid 
at least over a short record of 20 to 30 minutes, the waves 
and wave induced particle velocities can be dealt with as 
stationary random phenomena, thereby allowing the utiliza- 
tion of standard statistical analysis procedures. 

B. REVIEW OP PREVIOUS WORKS 

The study of the kinematics of breaking waves in the 
surf zone has progressed slowly both due to the difficulty 
in making direct field measurements, and due to the diffi- 
culty in modeling the surf zone in the laboratory. 

Early work by Iversen (1953) involved use of a wave 
channel in which waves were generated using a hinged plane 
oscillating flap. Measurements ’were made using photographic 
techniques. The length of the channel limited the wave type 
to plunging and surging breakers. Adeyemo (1970) made simi- 
lar laboratory measurements using hydrogen bubble and photo- 
graphic methods. 

Inman (1956) made field measurements of water particle 
velocity by measuring the drag force under the wave in order 
to infer particle velocities. Miller and Zeigler (1964) used 
both acoustic and electromagnetic current meters to determine 
the particle motion in the surf zone. Walker ( 1969 ) made 
studies using propeller type flow meters. Thornton (1968) 
and Steer (1972) made use of electromagnetic current meters 
to measure the particle motion. 



9 



C. OBJECTIVE 



The objective of this research is to study the kinematics 
of water particles in breaking waves within the surf zone. 
Simultaneous measurements were made of the instantaneous 
water surface elevation and two orthogonal water particle 
velocities at a fixed location in the surf zone. Additional 
wave measurements were made at points seaward and shoreward 
of the primary location. Estimates of the probability 
density functions and power spectra of the wave height and 
particle velocities are to be made. The applicability of 
linear wave theory is to be measured in computing velocity 
spectral components using the energy-density spectrum of 
the waves. The computed velocity spectra are to be compared 
with actually measured velocity spectra. 
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II. STATISTICAL ANALYSIS 



It was assumed that the data to be analyzed represents 
realizations of stationary, random processes. This assump- 
tion allows the use of standard statistical techniques. The 
basic analyses are described in the following sections as 
if performed on continuous data of infinite duration for 
ease of discussion. The actual analyses were performed 
digitally on records of finite length. 

A. PROBABILITY DENSITY FUNCTION 

The probability density function describes the probability 
that a sample time history record x(t) will assume a value 
between x and (x + Ax), and that this approaches an exact 
probability as T, the observation time, approaches infinity, 
and as the defined range, Ax, approaches zero. The proba- 
bility density function is given mathematically as: 

T 

P(x) = Lim Lim ~ ~ (2.1) 

Ax+0 T-* 00 x 1 

Where 1\ is the amount of time that x(t) falls within the 
range (x, x + Ax). 

The computed probability density function was compared 
to the Gaussian and Gram-Charlier distributions and tested 
for comparison using the chi-square goodness-of-f it test. 
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The Gaussian probability density function is given as: 



P 



Ga 



(x) 



1 e -(x-y) 2 /2a 2 

a V2F 



where a = standard deviation 
y = mean 



( 2 . 2 ) 



and in standardized form, with Z = , 



P (z) = l e -Z /2 

Ga \^T 



The Gaussian distribution is completely described by the 
mean and variance, and has zero skewness and kurtosis equal 
to three. 

The Gram-Charlier distribution was shown by Longuet- 
Higgins (1963) to better describe the probability distribu- 
tion of the instantaneous water surface elevation, f|, in 
moderately non-linear swell in deep water. It is given by: 



P Gc (z > “ P Ga (z)[1 '° + 2* h 3 



m 



4 

2? 



m. 



+ (^rir Hi, + 



ir V + 



...] ( 2 . 3 ) 



where 



it _ n n(n-l) z n 2 , n(n-l) (n-2) (n-3) z n ^ 
n z 1! 2 2! 2 2 

(Hermite Polynomial) 

m^ = skewness of record 

mjj = kurtosis - 3 
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It can be seen that for zero skewness and kurtosis of three, 
the fourth order Grara-Charlier distribution degenerates to 
the Gaussian distribution. 

To apply the chi-square goodness-of-fit test, the random 
variable, x(t), was divided into K intervals, called class 
intervals, thereby forming a frequency histogram. The 
observed frequency, f ^ , is the number of observations within 
the ith class interval. The expected frequency, F^, is the 
number of observations which fall in the ith interval of the 
theoretical distribution. The total discrepancy between the 
observed and expected frequency is given by: 




(2.H) 



2 

As the value of x increases from zero, the goodness-of- 



fit decreases. 



B. SPECTRAL ANALYSIS 



1 . Energy-Density Spectrum 

The auto-covariance for a record, x(t), is defined 



mathematically by 



, t+T 

R, -, ( t ) = Lim =• / x n (t)x, ( t+x)dt 

xx 1 t X x 



( 2 . 5 ) 



The value obtained describes the dependence of the values of 
the data set at one point in time on the values at some 
later time* t, in the data set . The function is real-valued and 



even, and has a maximum at x = 0 . The above expression 
implies an infinitely long data set. To deal with the 
reality of a data set of finite length, a lag window is 
applied. The choice of lag window for this analysis was the 
Parzen lag window, because it has the advantage of no 
negative side lobes, which maintains the coherence function 
between its theoretical values of ± 1. It also insures 
no numerical resonance problems when computing cross-spectra. 
The Parzen window is given by Bendat and Piersol (1966) as: 



1 - + 6 (J-) 



P(t) - 



m 



m 



2 Cl - =?-]■ 

m 

0 



, T = 0,1,2, . . 

T 

_ m , , 

5 T “ 2 + ’ 



, x > T 
5 m 



T 

m 

5 2 






( 2 . 6 ) 



m 



The energy-density spectrum is computed from the Fourier 
transform of the covariance function as modified by the 
Parzen lag window. The spectrum is given by: 

S 11 (f) = / P(x)R i;L (x)e- i2TTfT dx (2.7) 

— 00 

where f is a particular frequency of interest. 

2 . Cross-Spectral Density 

The cross-covariance function describes the general 
dependence of the values of one set of data on the other set, 
in contrast to the auto-covariance function which compares a 
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record against itself. The cross-covariance function is 
given mathematically as: 

t+T 

R (x) = Lim - / x (t)x ? (t+x)dt (2.8) 

ld r p+oo 1 ± d 

t/ 

The quantity R^ 2 (x) is a real-valued function, but it is not 
an even function, and does not necessarily have a maximum 
value at x = 0. 

A lag window again is applied to the cross- 
covariance function before computing the cross-spectral 
density. Applying the Parzen window and taking the Fourier 
transform, the cross-spectral density is given by: 

oo 

S 12 (f) = / P(T)R 12 (T)e" 12TrfT dr (2.9) 

—OO 

Since, unlike the energy-density spectrum, the cross-spectrum 
is not an even function, it is expressed in real and 
imaginary parts by: 



S 12 (f) = C 12 (f) " 1 Q 12 (f) 



( 2 . 10 ) 



The real part, called the co-spectrum, is given by: 



C 12 (f) =2 / P(t)[R 12 (t) + R 12 (-T)]cos(27Tfx)dT (2.11) 



The imaginary part, called the quadrature spectrum, is given 



by: 



Ql 2 (f) =2 / P(x)[R 12 (x) + Rj_ 2 (~x) ]sin( 2xfx)dx (2.12) 
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3 . Phase Angle 

Equation (2.10) can be rewritten in complex polar 
notation in terms of the modulus of the cross-spectral 
density and a phase angle, 



where = avera S e angular difference by which the 

cross-correlated components of x 2 (t) lead those of x^(t) in 
each spectral band. The phase angle expressed in terms of 
the co-spectra and quadrature-spectra is given by: 



4 . Coherence 

An expression for the comparison of two random 
variables as a function of frequency is described by the 
coherency squared, or simply coherence, and is derived using 
the expressions previously stated for energy-density spectra 
and for cross-spectral density. In mathematical terms, the 
coherence function is given as: 



S 12 (f) = l S i2 (f) l e iei2(f) 



(2.13) 




(2.14) 



where 



-it £ e 12 (f) - 7r 




y 12 (f) = 

Sii(f)S 22 (f) 



(2.15) 



where : 
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The coherence function has the characteristics that when 
2 

Yl 2 (f) = 1, the two records are completely correlated at 

2 

that frequency; and when (f) = 0, the records are 

statistically independent at that frequency. 
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III. EXPERIMENTAL PROCEDURE 



A. EXPERIMENTAL SITE 

The experiment was conducted at Del Monte Beach inside 
Monterey Bay adjacent to Monterey, California, on 12 April, 
1973 • The beach profile shown in Figure (1) is one of gentle 
slope. Plunging and spilling breakers are very common at 
this site, with spilling breakers being the predominant type. 
The waves at this location are usually severely directionally 
filtered by refraction due to the geometry of the coastline 
and bay, and approach almost perpendicular to the shore. 

B. INSTRUMENTATION 

The primary measuring instruments were placed on a three 
meter steel tower as shown in Figure (2). The tower was 
guyed so that tower movement and vibration were negligible. 
The instruments were arranged so as to be in the same verti- 
cal plane, with the exception of a small offset parallel to 
the wave crests given to the water current meter to prevent 
interaction with the capacitance probe. 

The water particle velocities were measured with an 
Engineering Physics Company water current meter Model 6130. 
The sea surface elevation was measured with a capacitance 
type penetrating wave staff fabricated locally, and with 
an Interstate Electronics Model DP-200 portable wave recor- 
der with SDP-201 differential pressure sensor. 
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FIGURE 1. Beach Profile at Experimental Site 




20 



FIGURE 2. Schematic of Main Instrument Tower 



Additional instrumentation was placed seaward and shore- 
ward of the main tower. A pressure wave gauge was installed 
seaward of the main tower by mounting it on a weighted sled. 
The sled was towed into position by attaching a winch cable 
from R. V. Acania. Acania, in her anchored position, heaved 
around, pulling the sled into position. The mean water at 
the sled during the experiment was 2.69 meters. 

A Baylor Company Wave Staff System Model - 13528R was 
installed 12.9 meters shoreward on a steel tower. This wave 
staff was placed in a line with the primary tower and the 
sled so as to be on an orthogonal with the predicted wave 
crests . 

The capacitance wave staff was designed by L. F. 

McGoldrick (1969) at the University of Chicago as an inex- 
pensive, but ultra-linear highly sensitive instrument for 
the detection of capillary waves in the 20 to 30 HZ range. 

It operates on the principal that changing the plate dimen- 
sions of a capacitor changes its capacitance and hence its 
influence on an electrical circuit. The probe was constructed 
of one-eighth inch 0. D. bronze wire; a dialectric of one- 
sixteenth inch wall thickness polystyrene tubing which 
fitted tightly over the w ire; and appropriate electrical 
connections. The sea water in which the probe was immersed 
constituted the second plate of the capacitor. The probe 
was held stationary on the tower, and as the sea level fluc- 
tuated, the capacitance changed accordingly. These fluctua- 
tions were sensed by a transistorized circuit operated by 
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two self-contained ten volt batteries. The electronics 
section was housed in a watertight case which was mounted 
on the tower. This enabled the connecting leads to be 
relatively short, thereby minimizing wire-to-wire capacitance. 
A static calibration of the system was made by immersing 
the staff in a water tank. The resulting calibration is 
shown in Figure (3)> 

The Interstate Electronics Pressure Sensor used on the 
primary tower and at the sled site is a small, unbonded 
strain gauge bridge. The sea pressure is coupled by a neo- 
prene diaphragm to a silicone fluid filling the interior. 

The arrangement is such that slow changes, such as tides, 
are lost through the hydraulic filter, and the transducer 
only senses the rapid pressure fluctuations of waves. The 
high-pass filter response time constant is 48 seconds. 

The Engineering Physics Company electromagnetic current 
meter measures two orthogonal components of water particle 
velocity through a range of zero to one and one-half M/SEC 
with a maximum output error of one percent of full scale 
reading. The instrument has two electrical time constant 
settings, 0.3 and 1.0 seconds. The operation of the unit 
is based on Faraday's principle of magnetic induction. 

Faraday discovered that currents are induced in a conductor 
which is put in motion through a magnetic field. This 
induction can be expressed by the vector equation: 



E = U X B 
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where e is the induced electric field, U is the velocity of 
->* 

motion, and B is the intensity of the magnetic induction. 

In the sea the induced electric field, E, is linearly pro- 
portional to the amount of fluid flowing through the magnetic 
field. The magnetic field, B, is most intense near the probe, 
and decreases with the distance from the probe according 
to the inverse square law. Thus, the flow velocity near the 
probe is weighted more heavily than that away from the probe, 
with only that water which is within a distance of two to 
three probe radii contributing significantly to the flow 
induced voltage. The radius of the probe is 3/8 inch. 

The flow meter was recalibrated using the method of Kra- 
pohl (1972) and Steer (1972). This procedure involves oscil- 
lating the probe in a water tank. The equipment arrangement 
is shown in Figure (4). The flow meter probe was mounted 
on a carriage which travelled back and forth on rails. The 
carriage was driven by a variable speed motor geared to an 
eccentric throw arm. The peak carriage velocity was calcu- 
lated to be the tangential velocity of the throw arm. The 
values of maximum carriage velocity, V' , and maximum measured 
velocity, V , were determined for different angular frequen- 
cies. Each pair of electrodes was calibrated by orienting 
them parallel to the flow. 

Figure (5) is a plot of the ratio of measured maximum 
velocity and maximum carriage velocity vs. frequency. The 
instrument's response is shown to be flat out to frequencies 
of about 0.5 Hz. At higher frequencies, the cart became 
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FIGURE *J. Electromagnetic Current Meter Calibration Assembly 
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unstable. In addition, waves were generated in the tank 
which affected the flow pattern. This upper frequency limi- 
tation was not deemed a significant factor in measurement 
of ocean waves in which the significant energy was found to 
be below 0.5 Hz. 

The Model 13529R Baylor Company Wave Staff System is an 
instantaneous water level measuring device consisting of two 
tensional half inch, 6x19 IWRC SS wire ropes. Together with 
these lengths of stainless steel wire rope, the transducer 
produces an electrically linear direct current output that 
is proportional to the amount of wave staff above a short 
circuit produced by the water surface . The working length 
of the staff, through previous modifications, had been 
shortened from 50 to 5 feet. The results of the recalibra- 
tion are shown in Figure (6). 

C. DATA COLLECTION 

The experiment was performed on 12 April, 1973- The 
towers and sled had been installed at low tide two days pre- 
viously. Tests were conducted and final calibration checks 
were made on 11 and 12 April. Included in the calibration 
was the placing of a two volt D. C. reference voltage on 
the FM tape to be used during data processing. This insured 
that the signal gains could be checked after digitizing. 

During the afternoon of 12 April, swell from a storm off 
the Washington coast began to arrive, with a buildup to peak 
wave height coincident with high tide. During the experimer , 
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the significant wave height at breaking was about 1 meter. 

The waves were generally of the plunging type, and were 
breaking in the general area of the main tower. The breaker 
point was visually observed to move seaward and shoreward 
according to the incident wave height. 

All data were simultaneously recorded using a Sangamo 
Model 3500, 1^ channel FM magnetic tape recorder. Continuous 
monitoring of data was accomplished using strip chart 
recorders and voltmeters. 

D. DATA PRE-PROCESSING 

The data were transcribed using the CS-500 Analog Compu- 
ter to an eight channel Clevit Brush strip chart recorder 
for signal verification. It was digitized and recorded onto 
seven track tape using the XDS-9300 computer in conjunction 
with the CS-500 Analog Computer. The digitized records were 
transferred to nine track tape for analysis by the IBM 360 
computer. 
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IV. ANALYSIS OP DATA 



Figures (7) and (8) are typical records showing that 
the breaking waves were roughly triangular in shape and were 
very peaked. This corresponds to the general observation 
that plunging breakers at or near the breaker point have a 
steep, almost vertical front face, and a sloping profile 
toward their trailing edge. The records show small peaks 
between the primary wave peaks . It is thought that the 
small peaks are harmonics of the primary wave frequency. 

The capacitance gauge showed a more sharply peaked wave 
than did the pressure gauges. This was due to the nature of 
the instruments in that the capacitance gauge could sense 
more rapid changes in the wave slope. 

A. SAMPLING PROCEDURE 

The data were read at a rate of 3 samples per second, 
corresponding to a sample interval of 0.33 seconds. This 
resulted in a Nyquist frequency of 1.5 Hz. The sampling 
rate was considered to be sufficiently high to avoid aliasing 
of energy into the portion of the spectra which was of 
interest from observations of the analog record. 

The length of record analyzed was 21 minutes. The maxi- 
mum lag time was taken as five percent of the record giving 
a spectral bandwidth resolution of 0.00873 Hz. According to 
Blackman and Tukey (1958), this yields 40 degrees of freedom. 
By applying the chi-square distribution, the 80 percent 
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FIGURE 7. Strip Chart Record of Waves 
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FIGURE 8. Strip Chart Record of Waves and Particle 
Velocities 
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confidence limits for ^0 degrees of freedom are found to 
be between 0.77 and I .38 of the measured power spectral 
estimates . 

B. PROBABILITY DENSITY FUNCTION 

The probability density function, pdf, for each wave 
and water particle velocity record was estimated using the 
procedures outlined in Section II-A. For each record, the 
standard deviation, variance, skewness, and kurtosis were 
computed and are summarized in Table I. Measured distribu- 
tions we re compared with the Gaussian and Gram-Charlier 
distributions using the chi-square goodness-of-f it test. 

The various measured distributions were also plotted against 
these analytical distributions. 

1 . Sea Surface Elevation 

At the seaward station, the measured wave pdf shows 
some relation to the Gaussian distribution [Figure (9)]. 

The Gram-Charlier gives a better fit in terms of the chi- 
square goodness-of-fit test. This is to be expected in 
that the Gram-Charlier pdf is calculated using the assigned 
parameters from the data, and including moments higher than 
the second . . 

At the primary station, the pdf again follows the 
Gram-Charlier more closely than the Gaussian, but not to as 
great a degree as in the deeper water case [Figures (10) 
and (11)]. 



33 









G 




















cd 








































O 




















1 


oa 


-=r 


o 


CJ\ 


.=4 


o 






cm e 


• 


• 


• 


• 


• 


• 








X cd 


VO 


OJ 


CO 


OO 


OJ 


1 — 1 








G 


t-- 


VO 


j=r 


t>- 


o- 


^r 








o 




rH 








OA 




















LfA 




























CO 


CO 


O- 


OJ 


OJ 


^r 


CO 






OJ C/3 


• 


• 


• 


• 


• 


• 








X 3 


OJ 


-=r 


CO 


OJ 


VO 


OA 








cd 


CO 


VO 


cr\ 


C-- 




CO 








C5 


CO 






rH 


i — i 


OJ 




c/3 








































o 




00 
















H 




H 


oa 




OJ 


j^r 


OJ 


OA 




Eh 




CO 


t— 


oa 




m 


oo 


CO 




< 




O 


VO 


lA- 


oo 


i — i 


CO 


t"- 




Eh 




Eh 


lta 


o 


-=r 


oo 


LfA 


o 




g> 




& 


• 


• 


• 


• 


• 


• 




Ph 




E> 


CO 


CO 




OJ 


OJ 






2 




X 
















o 


• 


















o 


Ph 




















• 


00 
















2 


Q 


CO 


VO 


1 — 1 




-=T 


1 — I 


p- 




o 


• 


w 


CO 


m 


o 


oa 


OJ 


OA 




M 


PH 




^r 


OJ 


CO 


OJ 


o 


rH 




EH 




Is 


VO 


VO 


1 — 1 


CO 


CM 


^r 




o 




pa 


• 


• 


• 


• 


• 


• 










o 


o 


rH 


o 


o 


o 




G> 




CO 
















pH 


















H 


>H 




t> 




















pa 


1 — 1 


o 


m 


GA 


OJ o 


VO o 


W 


H 




Q 


o 


OJ 


-=3* 


OJ 


O 0 


OO 0 


v-Q 


00 






OJ 


OJ 


i>- 


t>- 


-^r c/3 


O CO 


PQ 


2 




• 


r — 1 E 


OJ E 


iH E 


OJ E 


LTA\ 


o \ 


< 


w 




Q 










• E 


• E 


Eh 


Q 




Eh 


o 


o 


o 


o 


o 


o 








00 
















>H 




















EH 




















H 




















G3 




pa 










OJ 


OJ 




M 




o 


■=r 


CO 




LfA 


OO 


oo ^ 




PQ 






j=r 


oa 


o 


-■=r 


rH O 


tA- o 




< 




< 


i — loj 


-=r oj 


COOJ 


C— OJ 


OA Q) 


O 0 




PQ 




M 


o g 


o E 


o E 


o E 


OJ 00 


O CO 




O 




& 




• 






• \ 


• \ 




P^ 




< 


o 


o 


o 


o 


° B> 


° vB 




Ph 




> 


• 














Pn 




















O 




















>H 














CO 


OJ 




p^ 














G\ o 


i — 1 o 




<C 










VO 


o 


O 0 


OO 0 








< 






oo 


OA 


OJ CO 


O CO 




§ 




pa 






• E 


• E 


• \ 


• \ 




G> 

oo 




2 






rH 


i — i 


O E 


o E 










O 




















G 


0 






g 


u 










3 


G 






0 


0 










CO 


G 


0 




-p 


4A 










CO CD 


CO 0 


a 0 


0 


rH 0 


0 










CD W 


CO bO 


G bO 


hO 


cd PZ 


s: 










U 3 


0 3 


cd 3 


G 


4-^ 


1 — 1 










(U (d 


G cd 


4^ Cd 


cd 


G p 


Cd 










C5 


pH O 


*H O 


G C5 


O G 


o G 










G 




O 


o 


tsl 0 


•rH 0 










CD 0 


G 0 


cd 0 


i — 1 0 


•rH Pi 


G 










O > 


•H \> 


CH > 


>> > 


g g 


U G 










3 cd 


oJ cd 


cd cd 


cd cd 


O 3 


0 G 










o IS 


J>J ^ 


o ;s 


PQ Is 


ig o 


> O 



3*i 



i 




35 



FIGURE 9- Frequency Distribution of Wave Pressure at Seaward Station 
Plotted with Gaussian and Gram-Cnarlier Distributions 
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FIGURE 10. Frequency Distribution of Wave Pressure at Main Tower 
Plotted with Gaussian and Gram-Charlier Distributions 
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FIGURE 11. Frequency Distribution of Wave Height Measured with 
Capacitance Wave Gauge Plotted with Gaussian and 
Gram-Charlier Distributions 



At the shoreward station, the distribution again 
more closely resembles the Gram-Charlier [Figure (12)]. 



The positive skewness measured at each station indi- 
cates a large number of small negative values and a higher 
probability of large positive values of wave height and 
velocity when compared to the Gaussian. This qualitative 
observation supports the notion that ' in shallow water, the 
wave crests are steep and narrow, and the wave troughs are 
flat and wide . 

2 . Water Particle Velocity 

The probability density functions shown in Figures 
(13) and (1*1) for water particle velocity again follow the 
Gram-Charlier distribution more closely than they do the 
Gaussian. The Gram-Charlier distribution was computed by 
carrying the series to include the fourth moments of the 
frequency distribution (Equation 2.3)* This appeared to 
be sufficient at the seaward station [Figure (9)]* However, 
due to the more non-linear situation in the surf zone, it 
appears that higher moments should be included. This would 
give a better fit and get rid of the humps in the tails of 
the distribution. 

C. SPECTRAL ANALYSIS 

The energy-density spectrum was calculated for each of 
the data records. Cross-spectra were then computed between 
variables, from which the coherence and phase were determined. 
Wave energy-density spectral components were computed from 
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FIGURE 13. Frequency Distribution of Horizontal Water Particle Velocity 
Plotted with Gaussian and Gram-Charlier Distributions 
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STANDARD DEVIATION 

FIGURE 14 Frequency Distribution of Vertical Water Particle Velocity 
Plotted with Gaussian and Gram-Charlier Distributions 



pressure energy-density components using linear wave theory. 
Linear wave theory was also used to transform the wave spec- 
tral components to spectral estimates of horizontal and ver- 
tical velocities. These calculated spectral values were 
then compared to measured particle velocity spectra. 

1 . Wave Height 

The capacitance wave gauge was taken as the primary 
wave height measuring instrument. The capacitance gauge gave 
a direct measurement of instantaneous sea surface elevation, 
and by taking the mean, a measure of the average sea surface 
elevation was obtained. 

The waves were also measured using a pressure wave 
gauge. To spectrally convert these pressure fluctuations to 
values of water surface elevation, a transfer function was 
applied. This function is based on linear theory, and is 
given as: 



= r 00 sh kh 

(wave) PS cosh « h+z > 



S(f) 



(4.1) 
(pressure) 



where: S(f) 

k 
h 
z 



= energy-density spectra 
= wave number 
= water depth 
= pressure sensor depth 



The energy-density spectra of both the capacitance 
and pressure wave gauges are shown in Figure (15)* The 
pressure spectrum has been converted to a wave spectrum 
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PHASE ANGIE COHERENCE LOG SPECTRAL ENERGY DENSITY (M^-SEC) 

-200 000 +200 0.0 1.0 - 5-0 - 3.0 - 1.0 



Main Pressure 
Wave Gauge 






0.30 



r 

0.45 



0.60 



0.75 



FREQUENCY (HZ) 



FIGURE 15. Spectra of Wave Height Comparing Main 
Pressure Wave Gauge and Capacitance 
Wave Gauge 
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using Equation 4.1. The spectra, which are characteristic 
of all the records, show a narrow banded energy-density peak 
at a frequency of 0.07 Hz., corresponding to a period of 14.3 
seconds. Also evident are peaks at 0 . 1 ') and 0.21 Hz. which 
appear to be harmonics of the primary peak at 0.07 Hz. These 
harmonics appear to be physical as was observed from the strip 
chart earlier. These harmonics probably also have energy 
contributions due to the Fourier computational technique. A 
low xrequency peak is also shown at about 0.02 Hz. corres- 
ponding to a period of 50 seconds. Energy at this frequency 
is commonly measured along this coast and is attributed to 
surf beat. 

A peak occurs in all spectra at 0.6 Hz. This is 
attributed to mechanical noise in the tape recorder. This 
peak is at a frequency outside the range of appreciable wave 
energy and does not affect the results of the experiment. 

The magnitude of the wave spectral components at 
the frequency of peak energy density calculated from the 
pressure record using linear wave theory was approximately 
60 percent of that measured by the capacitance wave gauge. 

The coherence values were very high, ranging to greater 
then 0.9 in the maximum energy portion of the wave spectrum. 
The phase spectrum shows the capacitance wave gauge measure- 
ments generally leading the pressure measurements by a few 
degrees, probably due to the ability of the capacitance 
gauge to better measure the curling portion of the wave. 
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The spectra of the offshore wave gauge and the capac- 
itance wave gauge were compared, as were the inshore and 
capacitance gauges [Figures (16) and (17)]. A dynamic pressure 
correction [Equation 4.1] was not applied to the seaward 
pressure gauge. The results show that the waves were 
shoaling throughout the process; that is, there was a gradu- 
al buildup of energy as the wave progressed from outside 
the breaker zone into the breaker zone . 

The greater variance at the inshore wave gauge com- 
pared with the wave gauge at the breaker line is attributed 
to energy reflection from the beach. The Baylor wave gauge 
tower was approximately 25 meters from the point at which 
mean water level intersects with the beach, and 12.5 meters 
shoreward of the mean breaker line where the flow meter and 
primary wave gauges were located. These distances correspond 
to one-half and one-quarter wave lengths, respectively, of 
a 13 second period progressive wave in a depth of 1.9 meters. 
This arrangement places a node at the breaker line and an 
anti-node at the Baylor wave gauge tower for a standing wave 
of this period. This 13 second period corresponds closely 
to the period of peak energy. 

The phase for a theoretical progressive wave is 
plotted with the measured phase for comparison in Figure (17). 
The phase relationship was calculated using linear wave 
theory. Deviations from the calculated progressive phase 
angle could be attributed to wave reflections. 
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PHASE ANGLE COHERENCE LOG SPECTRAL ENERGY DENSITY (M^SEi 

-200 000 +200 0.0 0.1 - 5.0 - 3.0 - 1.0 



■. Seaward 
Pressure 
V. T ave Gauge 




0.15 0.30 0.45 0.60 0.75 





FIGURE 16. Spectra of Wave Height Comparing Seaward 
Pressure Wave Gauge and Capacitance 
’Wave Gauge 
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FIGURE 17. Spectra of Wave Height Comparing Shoreward 
Resistance (Baylor) Wave Gauge and 
Capacitance Wave Gauge 
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2 . Sea Surface Elevation and Water Particle Velocity 



The spectra of sea surface elevation and viator 
particle velocities were compared by using linear wave 
theory to convert from vrave spectral components to water 
particle velocity spectral components. The water surface 
elevation described by linear theory is given by: 

p = a cos (kx-at) 0*2) 

whe re 

p = sea surface elevation 
a = wave amplitude 
k = wave number 

x = horizontal cartesian coordinate 
o = frequency 
t = time 

The velocity components are given by: 



Horizontal velocity - u 



aa 



cosh k(h -MO cos(kx-ot) (4.3) 



sinh kh 



Vertical velocity 



w = aa 



sinh k(h fjO sln(kx-at) (4.4) 
sinh kh 



The horizontal velocity spectrum calculated from the wave 
spectrum is given by: 



S(f) 



( vel) 



[0 



cosh k(h + z) -j 2 
sinh kh 



(wave) 



(4.5) 



and the vertical velocity spectrum by: 
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S(f) 



( vel) 



= [a 



sinh k(h + z) -.2 . 

sinh kh J 



( wave ) 



(4.6) 



where the terms in brackets are the corresponding transfer 
functions . 

a. Horizontal V/ater Particle Velocity 

Linear theory under-calculates the horizontal 
velocity spectrum by 50 percent at the peak energy point 
[See Figure ( 18 )]. The values of coherence are high, 
ranging above 0.75* The phase angle, which according to 
linear theory should equal zero, ranges up to ^13 degrees 
from its average value of about 20 degrees. This shows the 
crest of the unstable breaking wave leading the horizontal 
particle velocities in the body of the flow. This compares 
well with the laboratory findings by Adeyemo (1970) that 
the maximum shoreward horizontal velocity for monochromatic 
waves occurs 25 degrees after the passage of the wave crest. 

b. Vertical Water Particle Velocity 

The measured rms vertical velocity was an order 
of magnitude less than the measured rms horizontal velocity, 
and showed very little wave induced motion [Figure (19)]. 

The coherence reached a maximum value of 0.37 at a frequency 
of 0.15 Hz., which is above the frequency of maximum energy 
as measured for the wave and horizontal velocity spectra. 

A phase angle of 130° to 1^5° was computed, in contrast to 
the theoretical value of 90° between wave height and vertical- 



velocity . 
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FIGURE 18. Spectra of Measured and Calculated 
Horizontal Velocity 
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FIGURE 19. Spectra of Measured and Calculated 
Vertical Velocity 
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The lack of coherence between water surface 
elevation and vertical velocity could be due to a noisy 
electrical signal from the instrument, a conversion of 
wave-induced motion to turbulent motion, or both. 
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V. CONCLUSIONS 



The probability density functions for the various records 
compared better with the Gram-Charlier distribution than with 
the Gaussian distribution when tested using the chi-square 
goodness-of-fit test. This result verifies the observed 
asymmetrical shape of the waves in the breaker zone and 
points out the importance of including higher order moments 
in describing wave phenomena in very shallow water. 

The values of the horizontal velocity energy-density 
spectral components calculated from wave spectra using linear 
theory indicate a qualitative, but not a quantitative rela- 
tionship. In general, linear theory underpredicts the 
magnitude of the velocity spectra. 

The high coherence between waves and horizontal velocity 
shows that most of the motion in the body of the breaking 
wave is wave-induced. A similar conclusion could not be 
drawn when comparing waves and vertical velocity. 

The spectral difference between water surface elevation 
and horizontal velocity was near that predicted by linear 
theory. The waves lead the horizontal velocity in spectral 
phase on the average by 20 degrees, implying that the 
" curling” crest of the wave arrives prior to maximum water 
particle velocity. 

In summary, the water particle motion in the breaking 
waves seems to be very much wave induced in the body of the 
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flow. Linear theory under-predicts the magnitude 
velocity spectra when computing from wave spectra. 
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